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Neutral pion photoproduction on the nucleon in a chiral quark model 
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A chiral quark-model approach is adopted to study the yp —> n^p and yn —> if'n reactions. Good descriptions 
of the total and differential cross sections and single-polarization observables are obtained from the pion pro¬ 
duction threshold up to the second resonance region. It is found that (i) the u = 0 shell resonance A(1232)/’33, 
the n = 1 shell resonances A(1535)5 u and A(1520)Di3, and the u = 2 shell resonance A(1720)/’i3 play crucial 
roles in these two processes. They are responsible for the first, second and third bump structures in the cross 
sections, respectively, (ii) Furthermore, obvious evidences of A(1650)Sn and A(1620)53i are also found in the 
reactions. They notably affect the cross sections and the polarization observables from the second resonance 
region to the third resonance region, (iii) The u-channel background plays a crucial role in the reactions. It has 
strong interferences with the ^-channel resonances, (iv) The f-channel background seems to be needed in the 
reactions. Including the f-channel vector-meson exchange contribution, the descriptions in the energy region 
Ey = 600 ~ 900 MeV are improved significantly. The helicity amplitudes of the main resonances, A(1232)F’33, 
A(1535)S 11, A^( 1520 )£)i 3 , A( 1720 )F’i 3 , A( 1650)5 u, and A( 1620)5 31 , are extracted and compared with the re¬ 
sults from other groups. 

PACS numbers: 13.60.Le, 14.20.Gk, 12.39.Jh, 12.39.Fe 


I. INTRODUCTION 


Understanding of the baryon spectrum and searching for the 
missing nucleon resonances and new exotic states are hot top¬ 
ics in hadronic physics H. Photoproduction of mesons is an 
ideal tool for the study of nucleon and A(1232) spectroscopies 
in experiments i]. Neutral pion photoproduction reactions 
are of special interest because the neutral pions do not couple 
directly to photons so that nonresonant background contribu¬ 
tions are suppressed (i.e., no contact term contribution) IH]. In 
the past few years, great progress has been achieved in experi¬ 
ments studying of the yp tt^p reaction at JLab, CB-ELSA, 
MAMI, and GRAAL. These experimental groups have carried 
out precise measurements of the differential cross sections 
and single-polarization observables with a large solid angle 
coverage and a wide photon energy range Odl- fl^ . Recently 
they also have finished some measurements of the double¬ 
polarization observables idllMl. Furthermore, in recent 
years significant progress has been achieved in experiments 
measuring the yn —» n^n reaction as well. In 2009, some mea¬ 
surements of the beam asymmetries for the yn —> n^n process 
were obtained by the GRAAL experiment in the second and 
third resonances region d]. Very recently, the quasi-free dif¬ 
ferential and total cross sections in the second and third res¬ 
onances region for this reaction were also measured by the 
Crystal Ball/TAPS experiment at MAMI IH]. Thus, improve¬ 
ment of the experimental situations gives us a good opportu¬ 
nity to study the excitation spectroscopies of the nucleon and 
A(1232). 

Stimulated by these new measurements, many partial- 
wave analysis groups, such as Bn Ga 1 20l - l^ . SAID 11231-25]. 
MAID d, Kent d, Julich dllland ANL-Osaka d, 
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have updated their analysis in recent years. For the yp — > 
reaction, good descriptions of the data up to the second 
and third resonances region have been obtained by different 
groups. However, the explanations of the reaction data and the 
extracted resonance properties from the reaction still exhibit 
strong model dependencies. For example, the yp couplings 
for some well-established resonances, such as V(1535)5'ii, 
V(1650)5 11 and V(1520)Di3, extracted by various groups dif¬ 
fer rather notably from each other. For the yn —> n^n reaction, 
consistent predictions from different mproaches can only be 
obtained in the first resonance region 1^. Because of the lack 
of data, the predictions from different models in the second 
and third resonances region are very different. Fortunately, in 
this energy region some new measurements of the cross sec¬ 
tion for the yn —> n^n reaction at MAMI IH] were reported 
about one year ago. 

These new data for the yn —> n^n reaction not only pro¬ 
vide us a good opportunity to extract more knowledge of the 
neutron resonances, but also shed light on the puzzle of the 
narrow structure around IT = 1.68 GeV observed in the exci¬ 
tation function of tj production off quasi-free neutrons by sev¬ 
eral experimental groups ll3l1 - [^ . This narrow structure has 
been listed by the Particle Data Group (PDG) as a new nucleon 
resonance V(1685) iH. However, many controversial expla¬ 
nations about this narrow structure, such as the V( 1650)5 n 
and V(1710)Pii coupled-channel effects, interference effects 
between V(1650)5ii, V(1710)Pn and V(1720)Pi3, and ef¬ 
fects from strangeness threshold openings, can be found in 
the literature 11351 - 1^ . In our quark model study, we find 
that the narrow structure around W - 1.68 GeV can be ex¬ 
plained by the constructive interferences between V(1535)5 n 
and V( 1650)5 n 1^ . Our conclusion is consistent with the 
analysis from the BnGa group ll^ |4^ . It should be men¬ 
tioned that, the yn coupling for V( 1650)5 n extracted by us 
and BnGa group has a positive sign, which is opposite to that 
of PDG 113^ . Now, two questions arise naturally: (i) Can 
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some clues about the controversially discussed A^(1685) be 
found in the yn —> reaction? (ii) Are the properties of 

N(1535)S 11 and ^(IbSO)^ n extracted from the rjN channel 
consistent with those extracted from the n^N channel? To bet¬ 
ter understand these questions, a systematic analysis of the 
recent data for the neutral pion production off nucleons is ur¬ 
gently needed. 

In this work, we carry out a combined study of the yp —> 
and yn —» reactions in a chiral quark model. By 

systematically analyzing the new data for neutral pion pho¬ 
toproduction on the nucleons, we attempt to uncover some 
puzzles existing in the photoproduction reactions and obtain 
a better understanding of the excitation spectra of the nucleon 
and A(1232). It should be mentioned that there are interesting 
differences between yp —> n^p and yn —> ji^n. In the yp reac¬ 
tions, contributions from the nucleon resonances of represen- 
tatio n [70,"^ 81 will be suppressed by the Moorhouse selection 
rule In contrast, all the octet states can contribute to 

the yn reactions. In other words, more states will be present 
in the yn reactions. Therefore, by studying neutral pion pho¬ 
toproduction on nucleons, we expect that the role played by 
intermediate baryon resonances can be highlighted. 

In the chiral quark model, an effective chiral Lagrangian is 
introduced to account for the quark-pseudoscalar-meson cou¬ 
pling. Since the quark-meson coupling is invariant under the 
chiral transformation, some of the low-energy properties of 
QCD are retained. The chiral quark model has been well de¬ 
veloped and widely applied to meson photoproduction reac¬ 
tions ll^l43l - l^ . Recently, this model has been successfully 
extended to ttN and KN reactions as well ifssl - l^ . 

The paper is organized as follows. In Sec. [Ill a brief review 
of the chiral quark model approach is given. The numerical 
results are presented and discussed in Sec. [Till Finally, a sum¬ 
mary is given in Sec. |IVl 


decay constant. The a>y is the energy of the incoming photons. 
The \Ni), \Nj) and \Nf) stand for the initial, intermediate and 
final states, respectively, and their corresponding energies are 
Ej, Ej and E y, which are the eigenvalues of the non-relativistic 
constituent quark model Hamiltonian H ll59l - l^ . The s- and 
M-channel transition amplitudes have been worked out in the 
harmonic oscillator basis in Refs, lliil - ld^ . 

The f-channel contributions of vector meson exchange are 
included in this work. The effective Lagrangians for the vector 
meson exchange for the ynV and Vqq couplings are adopted 
as llH] 


£y.v = (5) 

llljl 

-Cvqq = gVqq^jiytl + (6) 

^triq 

where A and V denote the photon and vector-meson fields, 
respectively; n stands for the 7r-meson field; gv^y and gvqq are 
the coupling constants. The f-channel transition amplitude has 
been given in the harmonic oscillator basis in Refs. 

It should be remarked that the amplitudes in terms of the 
harmonic oscillator principle quantum number n are the sum 
of a set of SU(6) multiplets with the same n. To obtain the con¬ 
tributions of individual resonances, we need to separate out 
the single-resonance-excitation amplitudes within each prin¬ 
ciple number n in the i-channel. Taking into account the width 
effects of the resonances, the resonance transition amplitudes 
of the i-channel can be generally expressed as jd^ll 




2Mp 


s - H- iMuYu 




(7) 


II. THE MODEL 


In this section, we give a brief review of the chiral quark 
model. In this model, the s- and M-channel transition ampli¬ 
tudes are determined by irilld^l] 


Ms = Yj^Nf\H,„\Nj){Nj\- -i- -He\Nd, (1) 

Mu - YiNflH,- -^- iNj)(Nj\H,„\Ni), (2) 

Ei — co,„ — Ej 


where //,„ and Elg stand for the quark-pseudoscalar-meson and 
electromagnetic cou pling s at the tree level, respectively. They 
are described by i44l - l46ll 


^ yi (3) 

j Jm 

He = -J]ejylA'^(Kr), (4) 

j 


where ij/j represents the y'-th quark field in a hadron, is the 
field of the pseudoscalar-meson octet, and /„ is the meson’s 


where -\/i = £, + cOy is the total energy of the system, a is the 
harmonic oscillator strength, is the mass of the s-channel 
resonance with a width Tr, and Or is the separated opera¬ 
tors for individual resonances in the s-channel. In the Chew- 
Goldberger-Low-Nambu (CGLN) parameterization, the tran¬ 
sition amplitude can be written in a standard form 0: 


Or - 


iffa -e + f^ 


.^(cr-q)o--(kxe) 


+ifj 


lql|k| 


lql|k| 
+ if! 


,R (o-■ k)(q ■ 6) .,R(o--q)(q- e) 


( 8 ) 


where cr is the spin operator of the nucleon, e is the po¬ 
larization vector of the photon, and k and q are incoming 
photon and outgoing meson momenta, respectively. In the 
SU(6)®0(3) symmetry limit, we have extracted the CGLN 
amplitudes for the i-channel resonances in the n < 2 shell 
for the yp —> jfip and yn —> uPn processes, which have been 
listed in Tables |I] and [III respectively. Comparing the CGLN 
amplitudes of different resonances with each other, one can 
easily hnd which states are the main contributors to the reac¬ 
tions in the SU(6)®0(3) symmetry limit. 
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TABLE I: The CGLN amplitudes of i-channel resonances in the n <2 shell for the yp jfip process in the SU(6)i8>0(3) symmetry limit. We 


have defined A = -( 


+ Dlql, B = 


+ I^A, C = 

AQ- 


= fiiii 


+ 4a, d = 


+ ^A,x^ 

50- 


= Ikllql 


P'liz) = P'l'iz) = ^. The ca^, m„. and 


fl, ■ 5a^--’ -- 3 ff 2 ’ ■ - dz ’ -■ / W - fl .2 • 

Ef stand for the energies of the incoming photon, outgoing meson and final nucleon, respectively. The is the constituent uoi d quark mass. 
l/Pq is a factor defined by l/p^ = 2/m,. Pi(z) is the Legendre function with z = cos0. 


resonance 

[A6,“+‘iV3,«,f] 

/f 

f2^ 


f! 

N(938)Eh 

[56,2 8^ Q Q] 

0 

2, giUq 

0 

0 

A(1232)P33 

[56,M0,0,0] 


A 


0 

A(1535)5u 

[70,2 8, 1, 1] 


0 

0 

0 

A( 1620)5 31 

[70,2 10,1,1] 


0 

0 

0 

A(1520)Oi3 

[70,2 8,1,1] 

-^'27^(1 +2^)^^ 

+i^k44APXz) 

0 

-i^fcMAF''(z) 

A(1700)Z)33 

[70,2 10,1,1] 



0 

+i^A:MAF''(z) 

A(1440)Pu 

[56,2 8^ 2,0] 

0 

.jjJJlllpjLCx 

36x18 19'^ma 

0 

0 

A(1710)Pu 

[70,2 8,2,0] 

0 

.jnJl±k-Erx 

^'36x18 

0 

0 

A(1750)P31 

[70,2 10,2,0] 

0 

- 11 V2 2 j k ^ 
'36x18 19 

0 

0 

A(1720)Ei3 

[56,2 8^2,2] 

+ 2t^)DPXz)x 

-i^^k-E-Dx 

90 12'^2/»„ 

-i^^kDP'^(z)x 

0 

A(1900)Pi3 

[70,2 8,2,2] 

-iWT2k(l + ^„)DP'Xz)^ 

-i-^Elk-^Dx 

'90 12'^2mr/ 

-/^iffeOP''(z)x 

0 

A(1985?)P33 

[70,2 10,2,2] 

+i^T2ki^-^)DP’Xz)x 

Dx 

'90 12'^6m« 

+i^^kDP'^(z)x 

0 

A( 1920)^33 

[56,10,2,2] 

0 

-i-^^kz^Dx 

90 9 Inii, 

^' 90 9 

0 

A( 1600)^33 

[56,10,2,0] 

'^^'w~k~CP2^z)x 

+i^^k-k-Cx 

-i^Tk^CP'iiz)^ 

0 

A( 1905)^35 

[56,^ 10,2,2] 




■2V2 31: /ID'YaF 

' 3 630„V, 

A(?)F35 

[70,2 10,2,2] 


+'Ti^^^3(z)23 

-i^AP''iz)x^ 

+'Ti^^3(z)^" 

A(1680)Fi5 

[56,2 8^ 2,2] 

+'■^(1 + 4;;)ap'Xz)2^ 


+i^-^AP'’{z)x^ 

-i^AP''(z)x2 

iV(?)Fi5 

[70,2 8,2,2] 

+'■^(1 + 4-)AP2iz)x^ 




A( 1950)^37 

[56,^^ 10,2,2] 

+'^Tt-AP’Xz)2^ 


— i2dl—!S—Ap"(y\Y^ 

^ 3 210rti„ 



Finally, the differential cross section dcrldQ., photon beam 
asymmetry S, polarization of recoil protons P, and target 
asymmetty T are given by the following standard expres¬ 
sions illillhl: 

dcr aea„{Ei + M^XE/ + Mff) 1 |q| -^ 2 

St '-• <«> 

4 

S = 2Re(//;//i -2 I///P, ( 10 ) 

/=1 

4 

P = -2lm(H;H2 + H;Hi)/Y,\Hi\\ (11) 

/=1 
4 

= 2lm(H;Hi+H;H2)/Yj\Hi\\ ( 12 ) 

1=1 


where the transition amplitudes //, in the helicity space can be 
expressed by the CGLN amplitudes fi llh^ : 


Hi 

= -=sm 0 cos-(/ 3 -H/ 4 ), 

V 2 2 

(13) 

H2 

= V2cos^[(/2-/i)H-sin2^(/3-/4)], 

(14) 

H3 

= ^sin 6 »sin^(/ 3 -/ 4 ), 

V 2 2 

(15) 

H4 

= V 2 sin ^ [(/2 H- /i ) H- cos^ ^ih- f a)]- 

(16) 

In Eq.(|9l), the fine-structure constant is well determined, 
and the ittAW coupling constantis related to the axial vector 
coupling gA by the generalized Goldberg-Treiman relation 


1 / gAMff _ 8„nn 

“ 47r i A j “ 47r ■ 

(17) 


T 
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TABLE II: The CGLN amplitudes of i-’Channel resonances in the n<2 shell for the yn —> T^n process in the 811(6)00(3) symmetry limit. 


resonance 


/f 

f? 

/f 

ff 

A(940)Pii 

[56,2 g Q Q] 

0 


0 

0 

A(1232)P33 

[56,M0,0,0] 




0 

A(1535)Sn 

[70,2 8,1,1] 

-'f *(1 + 

0 

0 

0 

yV(1650)S,i 

[70,8,1,1] 


0 

0 

0 

A( 1620)531 

[70,2 10,1,1] 


0 

0 

0 

A(1520)Oi3 

[70,2 8,1,1] 



0 

-i^AxP'^iz) 

A(1700)Oi3 

[70,^ 8,1,1] 


+'^5^'4xP'(z) 

0 


A( 1700)033 

[70,2 10,1,1] 



0 

-H^AxP'Tz) 

A(1675)Oi5 

[70,^ 8,1,1] 



-'-r2^^2:z 

■'■'■6"30;7i"^2^^2d) 

A(1440)Pn 

[56,2 g^ 2,0] 

0 

, ■ 47 V2 10. k 
^' 36x108 

0 

0 

A(1710)Pn 

[70,2 8,2,0] 

0 

, • 47V2 2 ,k f-. 

36x108 

0 

0 

A(1750)P31 

[70,2 10,2,0] 

0 

i ^^252 1 r, k (-• 

' 36x108 11 

0 

0 

NC)Pn 

[70,8,2,2] 

0 

• 47V2 Ik k n 
' 36x108 9 

0 

0 

A( 1910)^31 

[56,10,2,2] 

0 

i ^^2 72 Dx 

' 36x108 

0 

0 

A(1720)P,3 

[56,2 g_ 2,2] 

“' T ^ ^ ® ^2 

108 2 tnui 

0 

0 

A(1900)P,3 

[70,2 8,2,2] 

+ ^)DP',(z)x 

-i^k-^Dx 

108 bnid 

-i^kDP'^(z)x 

0 

(V(?)P,3 

[70,8,2,2] 

0 

-i^^k-^Dx 

' 108 2 6/n« 

2^^^^2 d)-*: 

0 

A(1985?)P33 

[70,2 10,2,2] 


-i^ik-^nx 

' 108 2 6m,1 


0 

A( 1920)^33 

[56,^ 10,2,2] 

0 

— Ox 

' 108 2 6mc, 


0 

A( 1600)^33 

[56,^^ 10,2,0] 


+i-^^k-^Cx 

108 2 6mci 


0 

(V(?)P,3 

[70,^ 8,2,0] 


+iJllk-^Cx 

M08 2^6mo^^2 

0 

yV(1680)Fi5 

[56,2 g^ 2,2] 


+'tI ^^^3^)23 

0 

0 

iV(?)Fi5 

[70,2 8,2,2] 


+'Tf5^^^3d)x2 

+ill5AP''(z)x2 

-i^5AP3(z)x2 

iV(?)Fi5 

[70,^ 8,2,2] 



TO d)d 

“'toT to 6^"^^3 d)x2 

A(?)F35 

[70,2 10,2,2] 





A(1905)T35 

[56,^ 10,2,2] 



, 12 A p//(' \ 2 

_;A2 12 _^4 p'Y73y2 

' 18 35 6m„ZlX3U;X 

A( 1950)^37 

[56,10,2,2] 


+'^5T(i7^^3d)v2 

;3 V2 8 A P^Y7') ^ 

' 4 9 210/»„ 


yV(1990)Fi7 

[70,^ 8,2,2] 



/3 V2 1 4 

^ 4 9 210/»„ ^-^4 



However, the quark model predicts rather large values gA - 
5/3 for charged pious and gA - 5 V2/6 for neutral pions. In 
our paper, the coupling a-„ is determined by htting the data. 

III. CALCULATION AND ANALYSIS 
A. Parameters 

In our framework, the s-channel resonance transition am¬ 
plitude, Or, is derived in the SU(6)®0(3) symmetry limit. In 
reality, the SU(6)®0(3) symmetry is generally broken due to, 
e.g., spin-dependent forces in the quark-quark interaction. As 
a consequence, conhguration mixings would occur. The con- 
hguration mixings break the SU(6)®0(3) symmetry, which 
can change our theoretical predictions. Furthermore, the he- 
licity couplings and strong decay couplings of some reso¬ 


nances might be over/under-estimated with the simple quark 
model. To accommodate the uncertainties in the symmet¬ 
ric quark model framework, we introduce a set of coupling 
strength paramete rs, Cr, for each resonance amplitude by an 
empirical method 115014^ : 

Or CrOr, (18) 

where Cr can be determined by htting the experimental ob¬ 
servables. In the SU(6)®0(3) symmetry limit one hnds Cr — 
1, while deviations of Cr from unity imply the SU(6)®0(3) 
symmetry breaking. 

In our previous study of the tj photoproduction on the nucle¬ 
ons, we found the conhguration mixings seem to be inevitable 
for the low-lying S -wave nucleon resonances A( 1535)5 n and 
A(1650)5ii, and D-wave nucleon resonances A(1520)Di3 
and A(1700 )Z)i 3. By including conhguration mixing effects 
in the S - and D-wave states, we explicitly express their tran- 
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sition amplitudes as follows; 

(19) 

The coefficients and can be related to the mix¬ 

ing angles. We adopt the same mixing scheme as in our pre¬ 
vious work 1^ . 


'5ii(1535)' 

/ 

,5n(1650). 

\ 


cos 0 s 
sin 6 s 


- sin 6 s 
cos 0s 


|70,2 8,1/2-) 
170,4 8 , 1 / 2 -) 


, (20) 


and 


' £>13(1520)' 


, £>13(1700), 

V 


cos 6 d 
sin 0 £) 


- sin 0 o 
cos &£) 


|70,2 8,3/2-) 
|70,4 8,3/2-) 


■ ( 21 ) 


Then, the coefficients defined in Eg.lfTOll are given by 


^s^u(?535) ^ cos0s(cos0s - sin6»s/2), (22) 

C^™’("i535) = /?4 sin 05 (sin 0s - 2 cos 0s), (23) 

*^s™a650) ^ sin 0s (sin 0s-H cos 0s/2), (24) 

*^s™a650) ^ cos 0s (cos 0s-H 2 sin 0s), (25) 

^ COS0Z)(COS0Z)-j-^sin0o), (26) 

<^E 520 ) = <sin0o(sin0B-2VIOcos0o), (27) 
‘-Oi°’(r700) ^ R 2 sm6D(sin6D + j^i^cos6d), (28) 
^Dntnoo) ^ ^4 Sd{cos 6d +2 VlO sin 6d). (29) 


The parameters R 2 and R 4 are introduced to adjust the over¬ 
all strength of the partial wave amplitudes of [70,^ 8] and 
[70,4 8], respectively, which may be overestimated or under¬ 
estimated in the naive quark model If f? = 1, one 

finds that the Cr parameters of 5 - and D-wave states can be 
explained with configuration mixings only. In the calcula¬ 
tion, the mixing angle between A^(1535)5' n and A^(1650)5 n 
is adopted to be 0s = 26°, determined in our previous 
work 11^ . Notice that in our work, we have adopted Isgur’s 
later conventions ll6^ where wave functions are in line with 
the SU(3) conventions of de Swart In this frame, we 
obtain a positive mixing angle 0s. However, in line with the 
old conventions of the SU(3) wave functions from Isgur and 
Karl’s early wor ks 115^ roOT . one obtains a negative mixing 
angle 0s IlSOl -l^ 1^ 1^ . This question has been clarified 
in Refs. l6^ . Furthermore, the mixing angle between 
N(1520 )Di 3 and N(1700)Di3 is adopted to be 6n ^ 10° as 
widely suggested in the literature l5lM53Ll^l^l^ . 

For the yp —> reaction, we obtain R^ - 1.0 and 

R 2 - 1.5 by fitting the 450 data points of differential cross 
section, and the 53 data points of total cross section col¬ 
lected in Tab. [Till While for the yn —> jiPn reaction, we ob¬ 
tain 0.7, Rf ^ 1.4 and ^ 1.0 by fitting the 

36 data points of total cross section around the second reso¬ 
nance energy region 1.30 < W < 1.72 GeV recently mea¬ 
sured at MAMI iH]. With these determined R parameters. 


from Fqs. l22ll29l one can obtain the overall strength param¬ 
eters and for the S - and D-wave resonances. 

The determined Cr values for these low-lying resonances 
are listed in Tab. m From the table, we find that to re¬ 
produce the data we need to introduce two large coupling 
strength parameters €^ 33 ( 1232 ) - 1-83 and € 533 ( 1620 ) - 1-8 
for A( 1232)^33 and A(1620)53i, respectively. The reason 
may be the well-known underestimation of their photocou¬ 
plings in the constituent quark model iiSlil. We also need 
to enhance the contributions of A(1520)Di3 by a factor of 
€z)i 3 ( 1520 ) - 1-4, which can not be explained with configu¬ 
ration mixings only. The underestimation of the resonance 
amplitude of A(1520)Di3 is also found in the yN rjN pro¬ 
cesses within the quark model framework 1^ . which is due 
to the underestimation of the photocoupling of A(1520)Z)i3 in 
the constituent quark model. In the tt® photoproduction pro¬ 
cesses, the data favor a smaller contribution of A(1535)5ii 
than that in the SU(6)®0(3) symmetry limit. In the yp rfip 
reaction, the strength parameter ^ 0.61 can be natu¬ 

rally explained with the configuration mixings between two S - 
wave nucleon resonances A(1535)5 n and A(1650)5 n with a 
mixing angle 0s - 26°. However, in the yn —> jiPn reaction, 
the strength parameter ^ 0.43 can not be well ex¬ 

plained with the configuration mixing effects only, we should 
introduce a parameter - 0-2 to slightly decrease the transi¬ 
tion amplitude of [70,^ 8 ,1/2-]. Furthermore, we find that in 
the yn jfin reaction, the enhancement of the contributions 
of A(1720 )Pi 3 might significantly improve the descriptions of 
the experimental data. 


TABLE III: 450 data points of differential cross section, and 53 data 
points of total cross section of yp —> included in our fits. The^^ 

datum point is aboul INdata =4.3. 


Data Refs. 

Obser. 

Ey 

Ndata 

T- 

Xi l^data 

[4] MAMI 

dcrldSl 

240,260,278 

27 

357 

13.2 

MAMI 

dcrldSl 

300 ~ 400 

112 

540 

4.8 

[9] CB-ELSA 

dcrIdQ. 

438 ~ 862 

311 

1204 

3.9 

[22 MAMI 

cr 

240 ~ 335 

20 

10 

0.5 

[9] CB-ELSA 

cr 

342 ~ 1138 

33 

37 

1.1 


To take into account relativistic effects, the commonly ap¬ 
plied Lorentz boost factor is introduced in the resonance am¬ 
plitude for the spatial integrals ll44|] . which is 

(7fi(k,q)->€3j(nk,r,q), (30) 

where y*. = M^v/F; and y^ - Mj^lEf. 

The nNN coupling a„ and the coupling goj„y ■ gojqq from oj- 
meson exchange in the t channel are considered as free param¬ 
eters in the present calculations. By fitting the experimental 
data of yp —> 71 ^p reaction (see Tab. HIDi. we get gj^NN - 13.2 
(i.e., a„ = - 13-8) and gojny ■ gojqq - 1-37. The nNN 

coupling determined in this work is compatible with the value 
gnNN — 13.5 adopted in other literature According 

to the decay of w ^ Try, one obtains g^j,„y ^ 0.32 Then 
the (xiqq coupling extracted by us is g^qg ^ 4.28, which is 
consistent with the value gi^qq ^ 3 suggested in Ref. iIt^ . 
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TABLE IV: The strength parameters Cr determined by the experi¬ 
mental data. 


Cr parameter 

yp TT^p 

yn —> 

^[70,A] 

Sii(1535) 
^[70,‘‘8] 
Sll(l535) 
^[70,78] 
Sii(1650) 
^[70,‘‘8] 
'-'511(1650) 
^[70,7 8] 
^Di3(I520) 
^[70,38] 

^013(1520) 

^[70,7 8] 

'-'013(1700) 

p[70,38] 

''013(1700) 

C/>33(1232) 

C53i(1620) 

Co, 5(1675) 

C/>i3(1720) 

^•^^-0.04 

n oq+0.04 
^•'^^-0.06 

1 41+0.15 

'•^3-0.09 

0.09 

1 C'3+0.02 

i-O-^_0.04 

l-80!°i° 

1.00 

1.00 

0 43+0 09 
'*■^■*-0.09 

-0.42 

n 97-f0.06 
-0.06 

1.12 

1 09-1-0.09 
-0.09 

-1.05 

0.08 

2.05 

1.83 

1.80 

1.00 

3.20)“;; 


There are another two parameters, the constituent quark 
mass niq and the harmonic oscillator strength a, from the tran¬ 
sition amplitudes. In the calculation we adopt their standard 
values in the the quark model, - 330 MeV and - 0.16 
GeV“. 

In the calculations, the n = 3 shell resonances are treated 
as degeneration; their degenerate mass and width are taken as 
M - 2080 MeV and T = 200 MeV, since in the low energy re¬ 
gion the contributions from the ;t = 3 shell are not significant. 
In the u channel, the intermediate states are the nucleon and 
A(1232) and their resonances. It is found that contributions 
from the ii > I shell are negligibly small and insensitive to 
the degenerate masses for these shells. In this work, we take 
Ml = 1650 MeV {M 2 = 1750 MeV) for the degenerate mass 
of n = 1 (n - 2) shell resonances. In the s channel, the masses 
and widths of the resonances are taken from the PDG or 
the constituent quark model predictions if no experimen¬ 
tal data are available. For the main resonances, we allow their 
masses and widths to change around the values from PDG ll^ 
in order to better describe the data. The determined values are 
listed in Tab. lYl As a comparison, the resonance masses and 
widths of both pole and Breit-Wigner parametrizations from 
the PDG iH are listed in Tab. [V] as well. It is found that 
the resonance masses and widths extracted by us are in good 
agreement with the values of pole parametrization. The rea¬ 
son is that when we fit the data, a momentum independent 
width is used, which is similar to the pole parametrization. 
It should be pointed out that V(1720)Pi3 seems to be a narrow 
state with a width of 120 MeV in our model, which is about 
one half of the average value from the PDG 1^ . However, 
our result is in good agreement with that extracted from the 
TT p ^ K^A reaction by D. H. Saxon et al. iH. The strong 
decay properties of V(1720)Pi3 will be discussed in detail in 
our another work. 

To know some uncertainties of a main parameter (Cr, Mr, 
r r) we vary it around its central value until the predictions are 
inconsistent with the data within their uncertainties. The ob¬ 
tained uncertainties for the main parameters have been given 


in Tabs. lIVI and IV] 


TABLE V: The masses Mr (MeV) and widths T^ (MeV) for the 
i'-channel resonances in present work compared with the values of 
Breit-Wigner (BW) and pole parametrizations from PDG 14 01. 


resonance 

(Mr, r^loura 

{Mr, r^|)pole 

{Mr, T^IbW 

A(1232)P33 

V(1535)Si, 

V(1650)Si, 

A(1630)S3i 

A(1520)Di3 

A(1720)Pi3 

1210)^,98)2 
1515)^,115)1“ 
1660);“, 150);“ 
1600)[;, 135);; 
1518);, 105)^0 
1685);“, 120)fn 

1210)|,100); 
1510)^,170);; 
i655)[;,i35);; 
i6oo)|;, i30)|; 
1510);, 110);“ 

1675+15 250+^^® 

1232);, 117); 
i535)|;, 150);; 
i655)|;,i40);; 
1630);;, i40)i; 
i5i5);,ii5)|; 
1720^^^ 250^^^^ 


B. yp —» T^p 

The chiral quark model studies of yp n^p were carried 
out in Refs, about twenty years ago. During the 

past two decades, great progress has been achieved for pion 
photoproduction at JLab, CB-ELSA, MAMI, and GRAAL. 
The new data sets are more accuracy and have larger solid 
angle coverage and wider photon energy range. The improve¬ 
ment of the experimental situations gives us a good oppor¬ 
tunity to test our model and study the excitation spectra of 
nucleon and A(1232) at the same time. All the intermedi¬ 
ate states in the s channel classified in the quark model with 
n < 2 are listed in Tab. J] It should be pointed out that in 
this reaction the contributions from the nucleon excitations 
with the representation [70,^ 8] are forbidden by the Moor- 
house selection rule ll^ 4l . In the n - Q shell, both nucleon 
pole and A( 1232)^33 contribute to the reaction. Comparing 
their CGLN amplitudes listed in Tab. |I] we can obviously see 
that A( 1232)7*33 plays a dominant role for its larger ampli¬ 
tudes. In the n - \ shell, two 5-wave states V(1535)5ii 
and A(1620)53 i and two D-wave states V(1520)Di3 and 
A(1700)D33 contribute to the reaction. Considering configura¬ 
tion mixing effects, we find that V(1650)5 n and V(1700)Di3 
can also contribute to the reaction. Similarly, from Tab.Uwe 
can find that V(1535)5ii and V(1520)Di3 play a dominant 
role in the ;t = 1 shell S -wave and D-wave resonances, respec¬ 
tively. In the n -2 shell eight P-wave resonances and five F- 
wave resonances contribute to the reaction. Comparing their 
CGLN amplitudes we find that V(1720)Pi3 and V(1680)7'i5 
play a dominant role in the n - 2 shell P-wave resonances and 
P-wave resonances, respectively. 

In present work, we have calculated the differential cross 
sections, total cross section, beam asymmetry, target asym¬ 
metry, and polarization of recoil protons from pion production 
threshold up to the second resonance region for the yp —> ji^p 
reaction. The model parameters are determined by fitting the 
450 data points of differential cross section from MAMI l|j,|5l] 
and CB-ELSA in the beam energy region 240 MeV< Ey < 
862 MeV, and the 53 data points of total cross section from 
MAMI iF^ and CB-ELSA in the beam energy region 
240 < Ey < 1138 MeV (see Tab. HIDi. The;^^^ datum point 
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is dibovA INdata — 4.3. Our results are compared with the data in FigsfTIfTOl 

I 


20 - 

10 - 

0 - 

30 

15. 

0 - 

16- 

8 

0 - 

5- 

0 - 

4- 

0 - 

4. 

0 - 


1 1 1 1 1 1 1 

" 200,1121 

1 1 1 1 1 1 1 

220,1137 

1 1 1 1 1 1 1 

240,1154 

1 1 1 1 1 1 I 

260,1170 

1 1 1 1 1 1 1 

278,115^ ^ 

^ I ' I ' I ' 

■ 300,1201 

32o'l217' ' : 

34o', 1232 ' : 

36o' 1247 ' : 

380,'l262' ' ; 

-1-1-1-1-1-1-1 

■ 400,1277 

-1-1-1-1-1-1-1 

420,1292 

-1-1-1-1-1-1-1 

438,1305 

—1—1—1—1—1—1—1 

466,1325 

—1—1—1—1—1—1—1 

488,1340 

I-'-1-'-1-'-1-! 

513,1358 

■ 1 ' 1 ' 1 ! 

538,1375 

' 1 ' 1 ' 1 ! 

560,1390 

-'-1-'-1-'-1-! 

587,1408 

611,1424 


■ 


/X- 


■ ' 1 ' 1 ' 1 ' 

639,1442 

1 1 1 1 1 1 1 

664,1458 

1 1 1 1 1 1 1 

688,1474 

1 1 1 1 1 1 T 

714,1474 

1 1 1 1 1 1 ^ 

738,1505 


790',1537' 

1 1 1 1 1 1 1 

815,1552 

1 1 1 1 1 1 1 

837,1566 

862,1580 

-1-1-1-1-1-1-1 

yv. 

-1-1-1-1-1-1-1 

—1—1—1—1—1—1—t 

—1—1—1—1—1—1—1 

-1-1-1-1-1-1- ^ 


w 


K 

A 

I 

I 

Q. 

a 

■C 

T3 


0 50 100 150 0 50 100 150 0 50 100 150 0 50 100 150 0 50 100 150 

(degree) 


FIG. 1: (Color online) Differential cross section of the yp —» n^p reaction as a function of scattering angle. Data are taken from 0] (solid 
circles), 0 (solid squares), (t^ (solid triangles), and [3] (open triangles). The first and second numbers in each figure correspond to the 
photon energy Ey (MeV) and the nN center-of-mass (c.m.) energy IV (MeV), respectively. 


The differential and total cross sections are shown in Figs.[T] 
and 13 respectively. It is seen that the chiral quark model can 
obtain a reasonable description of the data in a wide energy 
region Ey - 200 ~ 900 MeV. To clearly see the contribu¬ 
tions from different resonances, we also plot the energy de¬ 
pendent differential cross sections in Fig. [3 One can clearly 
see three bump structures in both the energy dependent dif¬ 


ferential cross sections and the total cross section. Accord¬ 
ing to our calculations, we find that A( 1232)^33 is respon¬ 
sible for first bump at Ey ^ 300 MeV. It governs the re¬ 
action in the first resonance region. Both N{1535)S n 
Af(1520)Di3 together dominate the resonance contributions 
in the second resonance region. They give approximately 
equal contributions to the second bump at Ey ^ 700 MeV. 
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W (MeV) 


FIG. 2: (Color online) Total cross section as a function of c.m. en¬ 
ergy W for the yp —> n^p reaction. Data are taken from 0] (open 
circles), 1?^ (solid circles), and (triangles). The results for 
switching off the contributions from A^(1650)Sii, A( 1620)5 31 , and 
t channel; and the partial cross sections for A(1232)f’33, A^(1535)5 u, 
A^(1520 )Di3, A^(1720)/’i 3 and the u channel are indicated explicitly 
by different legends in the figure. 

The Af(1720)F’i3 resonance might be responsible for the third 
bump at Ey ^ 1000 MeV. It should be mentioned that although 
A( 1620)5 31 and V(1650)5ii do not give obvious structures 
in the cross sections, they are crucial to give a correct shape 
of the differential cross sections from the second resonance 
region to the third resonance region (see Fig. |4]i. Switching 
off their contributions one can see that the total cross sections 
around Ey - 700 ~ 1000 MeV are overestimated slightly 
(see Figs. |2]i. The M-channel background plays a crucial role 
in the reaction, it has strong destructive interferences with 
A(1232)F33, V(1535)5 n and V(1720)Pi3. By including the 
f-channel vector-meson exchange contribution, we find that 
the descriptions of the cross sections in the energy region 


The beam asymmetries S in the energy region Ey = 220 ~ 


Ey - 600 ~ 900 MeV are improved notably, while without 
the f-channel contributions, the cross sections are underesti¬ 
mated obviously (see Figs.|2and|4l). Finally, it should be men¬ 
tioned that our quark model explanation of the first and second 
bump structures in the cross sections are consistent with that 
of the isobar model lUi] . However, our quark model explana- 



E (MeV) (MeV) 


FIG. 3: (Color online) Photon energy dependent differential cross 
section of the yp —> n^p reaction. Data are taken from il (left 
triangles), III (open squares), It^ (solid stars), It^ (open trian¬ 
gles), (open circles), (HI] (solid squares), HI] (solid circles), jH) 
(down triangles), Ol] (up triangles), (H (diamonds). The par¬ 
tial cross sections for A(1232)F’33, Af(1535)5n, V(1520)Di3 and 
Af(1720)/’i3 are indicated explicitly by different legends in the fig¬ 
ure. 


tion of the third bump structure differs from that of the isobar 
model ill]. In Ref. 01 , the authors predicted that the third 
bump might be due to three major contributions: A(1700)D33, 
^(1680)7^15 and V(I650)5 n, rather than V(1720)/’i3. Thus, 
to clarify the puzzle about the third bump structure in the cross 
section more studies of the reaction yp n^p are needed. 


900 MeV are shown in Fig. H In this energy region, the 
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®c.m. (degree) 


FIG. 4: (Color online) Effects of 7V(1650)Sn, A(1620)53i and t 
channel on the differential cross sections of t he y p —» 7 T°p pro¬ 
cess. Data are taken from H (open squares) and fTHl (solid squares). 
The results by switching off the contributions from A^(1650)5ii, 
A( 1620)5 31 , and t channel are indicated explicitly by different leg¬ 
ends in the figure. The first and second number in each figure cor¬ 
respond to the photon energy Ey (MeV) and the nN center-of-mass 
energy W (MeV), respectively. 


The polarizations of recoil protons P are shown in Fig. |7] 
In the low energy region Ey < 650 MeV, only a few old data 
with limited angle coverage were obtained. Recently, some 
precise new data in the higher energy region Ey ^ 700 ~ 
900 MeV were reported by the CBELSA/TAPS Collabora¬ 
tion From Fig. |7] it is found that our quark model 

descriptions are in reasonable agreement with the measure¬ 
ments in a fairly wide energy region Ey - 280 ~ 800 MeV. 
Above the photon energy Ey ^ 800 MeV, the quark model de¬ 
scriptions at both forward and backward angles become worse 
compared with the data. To clearly see contributions from 
different resonances, the energy dependent P at six angles 
6 'c.m. = 40°, 60°, 90°, 110°, 130°, 150° are shown in Fig. [8] as 


The target asymmetries T are shown in Fig. |9l Below 
the photon energy Ey ^ 700 MeV, only a few old data 


polarized data are not as abundant as those of differential 
cross sections. In the low energy region Ey < 600 MeV, 
until now no data on S at the forward and backward angles 
had been obtained. From Fig. |5] it is seen that the chiral 
quark model has achieved good descriptions of the measured 
beam asymmetries S in the energy region Ey = 220 ~ 800 
MeV. In the higher energy region Ey > 800 MeV, it is found 
that the chiral quark model poorly describes the data at for¬ 
ward angles. To clearly see contributions from different reso¬ 
nances, the energy dependent beam asymmetries at six angles 
dc.n,. = 20°, 60°, 90°, 125°, 150°, 170° are shown in Fig. i|as 
well. From the figure, it is found that the beam asymmetry S is 
sensitive to A(1232)/’33. Its strong effects not only exist in the 
first resonance region, but also extend to the second resonance 
region. If we switch off the contributions of A( 1232)^33, 
the beam asymmetry S changes drastically. Furthermore, we 
find that both V(1520)Di3 and V(1535)5n have strong ef¬ 
fects on the beam asymmetry Z around the second resonance 
region (i.e., Ey ^ 700 MeV), and without their contribu¬ 
tions, the beam asymmetry S in this energy region changes 
notably. In the higher energy region Ey > 800 MeV, it is found 
that the resonances A(1232)f’33, V(1520)Di3, V(1535)5'ii, 
V(1650)5'ii, A(1620)5'3 i and V(1720)Pi3 together with the 
M-channel background have equally important contributions 
to the beam asymmetry S. It should be mentioned that when 
the beam energy Ey > 800 MeV, many P- and F-wave states 
in the n - 2 shell begin to have obvious effects on the beam 
asymmetry Z as well. Thus, so many equal contributors in this 
higher energy region make difficult descriptions of the beam 
asymmetry S. 


well. It is found that an obvious dip structure appears around 
Ey - 700 MeV, which can be well described in the chiral 
quark model. The dip structure is due to the strong effects of 
A(1232)F33. When we switch off its contribution, we find that 
the dip structure disappears. Furthermore, from Fig. |8]it is 
obviously seen that the polarization of recoil protons P is sen¬ 
sitive to V(1520 )Di 3 and V(1535)5' n around the second res¬ 
onance region (i.e., Ey ^ 700 MeV). In the higher energy re¬ 
gion > 800 MeV A(1232)P33, V(1520 )Di 3, V(1535)5ii, 
V( 1650)5 11 , the M-channel background and other higher par¬ 
tial waves have approximately equal contributions to P, which 
leads to a complicated description of the higher energy data. 


with a very small angle coverage were obtained. Recently, 
some precise data with larger angle coverage in the higher 
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FIG. 5: (Color online) Beam asymmetry of the yp —> ;r°p process as a function of scattering angle. The data are taken from @] (solid left 
triangles), (open squares), Isill (open up triangles), 1 ^ (solid circles), (H (solid squares), Il3l (diamonds), 1 ^ (solid down triangles), 
01 (stars). The first and second number in each figure correspond to the photon energy Ey (MeV) and the kN center-of-mass energy W 
(MeV), respectively. 


energy region Ey ^ 700 ~ 900 MeV were published by 
CBELSA/TAPS Collaboration ifl^ . By comparing with the 
new data we find that our chiral quark model calculation ob¬ 
viously underestimates the target asymmetry T in the higher 
energy region Ey > 700 MeV, but the predicted tendency is 
in rough agreement with the data. In the low energy region 
Ey ^ 280 ~ 450 MeV, the data can be well described in the 
chiral quark model, though the data at forward and backward 
angles are still absent. In the energy region Ey ^ 450 ~ 660 
MeV, our quark model results are obviously smaller than the 
data at the forward angle. To test our model, we expect that 
more precise measurements with large angle coverage can be 
carried out in the energy region Ey < 700 MeV in the fu¬ 
ture. To clearly see contributions from different resonances, 
the energy dependent target asymmetries T at six angles 
0c.,n. = 27°,50°,83°,100M20M45°areshowninFig.[in]as 
well. The data show that there is a dip structure at the an¬ 
gle 9c.m. - 80° ~ 100° around the second resonance region 
Ey - 700 MeV. This structure can be explained by the strong 
interferences between A(1232)P33 and V(1535)5n. Switch¬ 
ing off the contributions of either A(1232)P33 or V(1535)5 n. 


no obvious dip structure around Ey - 100 MeV can be found 
in the target asymmetry T. According to the chiral quark 
model predictions, the dip structure should be found at the 
forward and backward angles as well. In the higher energy 
region Ey > 700 MeV, it is found that many contributors, 
such as A(1232)P33, V(1535)5ii, V(1650)5ii, V(1520)Di3, 
A(1620)53 i, V(1720)Pi 3 and the w-channel background have 
obvious effects on the target asymmetry T. 

In brief, obvious roles of the A(1232)P33, V(1535)5ii, 
V(1650)5ii, A(1620)53 i, V(1520)Di 3 and V(1720)Pi3 have 
been found in the yp —» jfip process, (i) A(1232)P33 not only 
plays a dominant role around the first resonance region, its 
strong contributions also extend up to the third resonance re¬ 
gion, which can be obviously seen in the cross section, beam 
asymmetry, target asymmetry and polarization of recoil pro¬ 
tons. (ii) Both V(1520 )Z)i 3 and V(1535)5 n play a dominant 
role around the second resonance region. They are the main 
contributors of the second bump structure in the energy depen¬ 
dent differential cross section and total cross section. Their 
strong effects on the polarization observables can be seen ob¬ 
viously as well, (iii) V(1720)Pi3 might play a crucial role 
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in the third resonance region. It might be responsible for the 
third bump structure in the energy dependent differential cross 
section and total cross section. However, no dominant role 
of A^(1720 )Pi 3 is found in the polarization observables. It 
should be pointed out that the evidence of A^(1720)Pi3 around 
the third resonance region should be further confirmed due to 
our poor descriptions of the polarization observables in the 
higher energy region, (iv) A( 1620)5 '31 and A^( 1650)5 n are 
crucial to give the correct shape of the differential cross sec¬ 
tions in the second resonance region, although they do not 
contribute obvious structures in the cross sections, (v) Fur¬ 
thermore, the u- and f-channel backgrounds play crucial roles 
in the reaction as well. The u channel has a strong interfer¬ 
ence with the resonances, such as A(1232)P33, A^(1520)Di3 
and A^(1535)5ii. Including the f-channel vector-meson ex¬ 
change contribution, we find that the descriptions in the en¬ 
ergy region Ey - 600 ~ 900 MeV are improved obviously, 
(vi) No obvious contributions of the other resonances, such as 
N(1700 )Di 3, A(1700)Z)33 and N( 1680 )Fi 5 , are found in the 
jp —> 71^p process. 


300 600 900 300 600 900 

(MeV) (MeV) 



FIG. 6: (Color online) Photon energy dependent beam asymmetry 
of the yp —> n^p reaction. The data are taken from (right tri¬ 
angles), (op en up triangles), ( open diamonds), (solid 
diamonds), 11 ^ (down triangles), (open circles), 11 ^ (solid 
squares), HI (solid circles), and (open squares). The results 
by switching off the contributions from various partial waves are in¬ 
dicated explicitly by different legends in the figure. 


C. yn 

The chiral quark model studies of yn —> n^n were carried 
out in Refs. 04^ 1^ 1^ about twenty years ago. However, 
the model studies were limited in the first resonance region, 
because only a few scattered data were obtained from the 
old measurements in the early 1970s. Fortunately, obvious 
progress has been achieved in experiments in recent years. In 
2009, some measurements of the beam asymmetries for the 
yn n^n process were obtained by the GRAAL experiment 


in the second and third resonances region S. In this en¬ 
ergy region, recently the quasi-free differential and total cross 
sections for this reaction were also measured by the Crystal 
Ball/TAPS experiment at MAMI |[3l. Thus, these new mea¬ 
surements in the higher resonances region provide us a good 
opportunity to extend the chiral quark model to study these 
high-lying resonances. 

The contributors of the s-channel intermediate states clas¬ 
sified in the quark model with n < 2 have been listed in 
Tab. M In the n - 0 shell, the dominant contribution to 
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FIG. 7: (Color online) Polarization of recoil protons P of the yp —> ;r®p reaction as a function of scattering angle. Data are taken from 
(circles), (triangles), and ca (squares). The first and second number in each figure correspond to the photon energy Ey (MeV) and the 
jiN center-of-mass energy W (MeV), respectively. 


the reaction comes from the A(1232)F’33, which has much 
larger CGLN amplitudes than the nucleon pole. In the n - 

1 shell, three ^-wave states Ai(1535)5ii, Ai(1650)5'ii and 
A(1620)5'3 i[ 70,^ 10], and four D-wave states Ai(1520)Di3, 
Ai(1700)Di3, Ai(1675)Di5 and A(1700)D33 contribute to the 
reaction. By comparing their CGLN amplitudes listed in 
Tab. HD we find that N(1535)S n and Ai(1520)Di3 play dom¬ 
inant roles in these S - and D-wave resonances. In the n = 

2 shell, twelve P-wave resonance and seven P-wave reso¬ 
nances contribute to the reaction. Most of the P-wave and 
P-wave resonances in the n - 2 shell have comparable am¬ 
plitudes. Ai(1720)Pi3, Ai(1900)Pi3 and A(1600)P33 have rela¬ 
tively larger CGLN amplitudes in the n - 2 shell P-wave res¬ 
onances, while Ai(1680)Pi5 and A(1905)p35 have relatively 
bigger CGLN amplitudes among the n - 2 shell P-wave res¬ 
onances. 

In this work, we have carried out a chiral quark model study 
of the yn —» n^n reaction up to the second and third reso¬ 
nances region. In the SU(6)®0(3) symmetry limit, the pa¬ 
rameters from the u- and f-channel backgrounds and the A 
resonances A(1232)P33 and A( 1620)5 '31 for the n^n channel 
should be the same as those for the n^p channel, which have 


been well determined by the yp data. Thus, in the yn —> n^n 
reaction these parameters are taken to have the same values 
as in the yp n^^p process. The other strength parame¬ 
ters, Cr, for the main resonances Af(1535)5 n, Af(1650)5n, 
Af(1520)Di3 and Af(1720)Pi3 for the yn reaction can not be 
well constrained by the yp data for their different photocou¬ 
plings, thus, we determine them by fitting the 36 yn data 
points of total cross section around the second resonance en¬ 
ergy region 1.30 GeV< W < 1.72 GeV recently measured at 
MAMI iH]. The;t'^ datum point is about/Ndata - 2.8. Our 
results compared with the data have been shown in Figs HTIfTbl 
The differential cross sections compared with the data are 
shown in Fig.[TT] In the energy region what we consider, only 
a few data can be obtained. Fortunately, the abundant data 
for the yp —» n^p process help us well constrain some im¬ 
portant model parameters, as we pointed out above. From 
Fig. [TT] one can see that the data of the yn —> jiPn reac¬ 
tion are reasonably reproduced. To clearly see the contri¬ 
butions from different partial waves, we plot the energy de¬ 
pendent differential cross sections in Fig. [12] as well. Our 
results obviously show three bump structures in the forward 
angle region. It is found that the resonances A(1232)P33, 
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FIG. 8: (Color online) Photon energy dependent polarization of re¬ 
coil protons for the yp —» T^p reaction. Data are taken from da 
(open circles), m (solid circles), (diamonds), (solid down 
triangles), m (solid stars), ll^ (open stars), 1 ^ (solid squares), 
(solid up triangles), (solid left triangles), 19^1 (open down 
triangles), and (solid right triangles). The results hy switching 
off the contrihutions from various partial waves are indicated explic¬ 
itly by different legends in the figure. 


Af(1535)5' 11 , Af(1520)Di3 and Af(1720)Pi3 play crucial roles 
in the yn n^n reaction. The A( 1232)7*33 resonance is re¬ 
sponsible for the first bump structure around Ey ^ 300 MeV. 
Both N{152)5)S n and Af(1520)Di3 are the main contributors 
to the second bump around Ey ^ 700 MeV. The V(1720)Pi3 
resonance is most likely responsible for the third bump around 
Ey ^ 1000 MeV. 

The total cross sections compared with the data are shown 
in Fig. [13 Obvious roles of A(1232)P33, V(1535)5'ii and 
V(1720 )Pi 3 in the yn —» n^n reaction can be found in 
the total cross section as well. Recently, the total cross 
section was measured by the Crystal Ball/TAPS experiment 
at MAMI H. There are two obvious bump structures in 
the cross section in the second and third resonances re¬ 
gion (see Fig. [T3l l. The bump structure around the second 
resonance region receives approximately equal contributions 
from V(1535)5' n and V(1520)Di3, while the bump structure 
around the third resonance region might be due to the con¬ 
tributions of V(1720 )Pi 3. There are no measurements of the 


total cross section in the first resonance region. In this energy 
region, we predict that the ratio of total cross section between 
the TiPn channel and the jf'p channel (Tnlo-p is around 1 (see 
Fig.m. 

Furthermore, by analyzing the data of differential and total 
cross sections, we find that A(1620)53i and V(1650)5 n play 
obvious roles around their mass threshold. If we switch off 
them, the cross sections around their mass threshold are over¬ 
estimated significantly. It should be mentioned that, the role of 
N{165Q)S 11 should be confirmed by more accurate data in the 
future, which will be further discussed in Sec. lIIIDl Finally, it 
should be pointed out that the backgrounds play a crucial role 
in the reaction. The i/-channel background has strong destruc¬ 
tive interferences with A(1232)7’33, N{1535)S n, V(1520)Z)i3 
and ^(1720)7*13. Including the f-channel vector-meson ex¬ 
change contribution, we find that the descriptions of the cross 
sections in the energy region Ey - 600 ~ 900 MeV are im¬ 
proved significantly. 

The polarization observations for the yn —> n^n reaction 
are very sparse. In the year of 2009, the beam asymmetry 2 
in the second and third resonances region was measured by 
the GRAAL Collaboration for the hrst time ifl^ . Our chi¬ 
ral quark model results are shown in Fig. [15] From the fig¬ 
ure, it is seen that the model results are in rough agreement 
with the data. Our results are notably smaller than the data at 
the intermediate angles. To clearly see the contributions from 
different partial waves, the energy dependent beam asymme¬ 
tries 2 at six angles Oc.m. - 20°, 52°, 91°, 123°, 144°, 163° are 
shown in Fig.[T 6 |as well. From the figure, one can find that 
below the photon energy Ey ^ 500 MeV, the beam asym¬ 
metry is sensitive to A(1232)7’33 and the u-channel back¬ 
ground. By turning off one of them, the beam asymmetry 
changes drastically in this energy region. Similarly, we can 
obviously hnd that around the second resonance region, i.e., 
Ey ~ 700 MeV, the V(1535)5ii, V(1650)5ii, V(1520)Di3, 
A( 1232 ) 7 ’ 33 , A( 1620 ) 53 i and the u-channel background have 
strong effects on the beam asymmetry. Up to the second 
resonance region, the higher partial wave states, such as 
V(1720 )Pi 3 begin to contribute to beam asymmetry. Many 
resonances together with the backgrounds have approximately 
equal contributions to the beam asymmetry, leading to a very 
complicated description of the data. 

As a whole, a reasonable chiral quark model description 
of the yn —> jf'n reaction is obtained from the pion produc¬ 
tion threshold up to the second resonance region. Obvious 
evidences of the A(1232)7’33, V(1535)5ii, V(1520)Di3 and 
V(1720 )Pi 3 are also found in the yn — > jf'n reaction, (i) 
The ground state A(1232)P33, the ^-wave state V(1535)5'n 
together with the D-wave state V(1520)Di3, and the P- 
wave state V(1720)Pi3 are responsible for the first, second, 
and third bump structures in the cross sections, respectively, 
(ii) Furthermore, another two 5-wave states A( 1620)5 31 and 
V( 1650)5 11 have obvious effects on the differential cross sec¬ 
tion around their mass threshold, although they do not give 
any structure in the cross sections. It should be pointed out 
that the role of V( 1650)5 n should be further confirmed in fu¬ 
ture experiments, (iii) The backgrounds play a crucial role in 
the reaction. The u channel background has a strong construe- 
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FIG. 9: (Colo r online) Target a symm etry of the yp —» ;r°p reaction as a function of scattering angle. The data are taken from ll^ (solid 
circles), llOOll (solid triangles), 1 10 ill (open circles), and fl^ (solid squares). The first and second number in each figure correspond to the 
photon energy Ej (MeV) and the nN center-of-mass energy W (MeV), respectively. 


live interference with the y-channel resonances A( 1232)^33, 
Al(1535)5'ii and Al(1520)Di3. By including the f-channel 
vector-meson exchange contribution, we find that the de¬ 
scriptions in the energy region Ey = 600 ~ 900 MeV are 
slightly improved, (vi) No obvious evidence of the other res¬ 
onances, such as V(1700 )Di3, V(1675)Di 5, A(1700)D33 and 
V( 1680 )Fi 5 , is found in the yn —» process. 


D. Helicity amplitudes 

The accurate data for the jn —> n^n and yp —> n^p pro¬ 
cesses provide us a good platform to extract the helicity am¬ 
plitudes of the dominant resonances in these reactions. The¬ 
oretically, the helicity amplitudes for a baryon resonance 
N* photoexcitation on a nucleon are defined by 

Ai = = A\He\N;J, = A-l}, (31) 

where T = 1/2 and 3/2. As we know, the helicity amplitudes 
of a resonance are related to the transition amplitudes of the 
photoproduction reactions. Thus, we can extract the helicity 
amplitudes from the neutral pion photoproduction processes 


by the relation 


At’.P 

^1/2,3/2 


IqlMjjTs 


(32) 


where = r„om/Tf> is the branching ratio of the resonance. 
The quantity ^ for different resonances can be analytically ex¬ 
pressed from their CGLN amplitudes. We have given the ex¬ 
pressions of the f for several low-lying nucleon and A res¬ 
onances in Tab. |Vll We estimate the helicity amplitudes for 
these main contributing resonances: A(1232)P33, A(1620)5' 31 , 
V(1535)5ii, V(1650)5ii, V(1520)Di3 and ^(1720)^13. The 
branching ratios for ^(1720)7*13 are adopted from our 
quark model prediction, and the branching ratios for other res¬ 
onances are taken from PDG14 ill (see table rvni i. Our ex¬ 
tracted helicity amplitudes are listed in Tab. IVIIII As a com¬ 
parison, in the same table we also show our previous solution 
extracted from the 77 photoproduction process es 11 ^ . the re¬ 
cent analysis of the yN data from SAID 123142511 . Kent 
and BnGa |^, the average values from PDG14 l34|] . 
and the theo retical predictions from different quark mod¬ 
els iM[i^. 


j 


j 
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FIG. 10: (Color online) Photon energy dependent targ et as ymme- 
try of the yp jfip reaction. The d ata are taken from llOOll (solid 
squares), 1 10 ill (open squares), Il02ll (solid circles), and 18^ (open 
triangles). The results by switching off the contributions from vari¬ 
ous partial waves are indicated explicitly by different legends in the 
figure. 


From Tab. IVIIII it is found that the helicity amplitudes of 
A(1232)f’33 extracted in present work are in good agreement 
with the values from PDG14J3 and other partial wave anal¬ 
ysis groups Ml2iMl^llO^. 

The and A"^2 of A^(1535)5 n extracted in this work are 
compatible with the PDG average values and the latest analy¬ 
sis of the jN data from the SAID ll23l - l^ . It should be pointed 
out that with the same model we found a smaller yp coupling 
A^^2 ^ 60 X 10-^GeV-‘/2 for the N(1535)S n by the analysis 
of the yp —>7^7 process 1 ^ . and similar solution was also 
obtained in ll 10|] . The reason for the different yp cou¬ 
plings for N{1535)S n in the n^p and rjp channels should be 
clarified in future studies. 

All the partial wave analysis groups have extracted simi¬ 
lar yp coupling A^^2 for A( 1650)5 n from the data, which is 
also con s isten t with the theoretical predictions in quark mod¬ 
els i[Iol[IH]. However, contradictory solutions for the yn 
coupling A"^2 of A(1650)511 are obtained by different groups. 
Our previous analysis of the yn —» rjii reaction indicates 
a positive helicity coupling - 24 x lO^^GeV^'^^ for 


A( 1650)511 which is supported by the latest analysis 

of the same reaction from the BnGa jlA and Kent 
groups. However, in present work by analyzing the recent the 
final-state-interaction (FSI) corrected data of the yn —> uPn re¬ 
action from the A2 Collaboration IH], a negative helicity cou¬ 
pling Aj ^2 - “18 X 10“^GeV“62 js obtained, which is com¬ 
patible with the values from the PDG 14 if^ and the recent 
SAID analysis ll23l - l^ . Contradictory results for the yn cou¬ 
pling Aj ^2 of A( 1650)511 obtained from two different reac¬ 
tions with the same model indicate that the A( 1650)5 n state 
found in the yn —> n^n is possibly not the same state found 
in the yn —> pn if the data are accurate enough. It should 
be noted that the FSI is rather rough correction that assumes 
identical effects on the proton and the neutron, which certainly 
does not have to be the case H. Thus, considering that the 
data from the A2 Collaboration might bear large uncertainties 
in the second resonance region, with a small positive helicity 
amplitude, A';^2 ^ 20 x lO^^GeV-62 for A( 1650)5 n, we pre¬ 
dict the differential and total cross sections around the second 
resonance region ( see Fig.fTTb. If A(I650)5 n has a positive 
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FIG. 11: (Col or online) Diffe rentia l cross sections of th e yn n^n reaction as a function of scattering angle. Data are taken from Q] (open 
circles), Il03ll (solid circles), ilOdll (solid squares), and llOSll (solid triangles). The first and second number in each figure correspond to the 
photon energy Ey (MeV) and the nN center-of-mass energy W (MeV), respectively. 


helicity amplitude, it is found that i) the differential cross sec¬ 
tion and the total cross section around the second resonance 
region should be significantly larger than the present data; and 
ii) Ai(1650)5'ii has obviously constructive interference with 
Ai(1535)5' 11 and Ai(1520)Di3, which can be tested in future 
experiments. It was pointed out in Ref. CH that the positive 
A "/2 would imply Ai( 1650)5 n should have a large ss compo¬ 
nent in its wave function. To clarify the sign problem of the 
yn coupling for A^( 1650)5 n, more accurate data are needed. 

We find a large helicity amplitude for A( 1620)5 31 , which is 
about a factor 2 larger than the PDG average valu^JM], and 
30% larger than the recent results from the BnGa 11^ 1^ and 
SAID groups. However, we find that our res ult is ver y 
close to the theoretical predictions in quark models IIIOTLIIOsIi . 


In our previous work 1^ . we gave our estimations of the 
helicity amplitudes for A(1520)Di3 by the analysis of the 77 
photoproduction data. However, the large uncertainties of the 
branching ratio lead to a weak conclusion of these helic¬ 
ity amplitudes. In this work, the accurate branching ratio b„N 
should let us extract the helicity amplitudes for A(1520)Di3 
more reliably. It is found that the extracted by us is in 
good agreement with the results from SAID group 12^ and 
the PDG average value 1^ . However, the A^^j extracted in 
present work are about 30% smaller than the PDG average 
value and the results from other groups. It should be men¬ 

tioned that recently the CBELSA/TAPS Collaboration also 
found a small helicity amplitude A^^j - 118 x 10 “^GeV“*^^ 
from an energy-independent multipole analysis based on new 
polarization data on photoproduction of neutral pions iI3. 
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FIG. 12: (Color online) Photon energy dependent differential cross 
sections of the yn —> Tfin rea ction . Data are taken from 1311 (so lid tri¬ 
angles), Il06ll (solid stars), Il04ll (solid squares), and TlOSll (solid 
circles). The partial cross sections for A(1232)f’33, A^(1535)5ii, 
N( 1520)Di 3 and N{\120)P\^ are indicated explicitly by different leg¬ 
ends in the figure. 


The yn couplings for the Af(1520)Di3 extracted in this work 
are compatible with the PDG values within 30% uncertainties. 
Our results are slightly smaller than the results from other par¬ 
tial wave analysis groups. 

For Af(1720)Pi3, we have noted that the absolute values of 
the ^nd extracted by us are compatible with the re¬ 
sults from the BnGa ll^l2^ and Kent 12711 groups. However, 
their solutions have opposite signs to our results. It is inter¬ 
esting to find that our results are consist ent w ith the quark 
model predictions by Z. Li and F. Close IIIOTII . and the par¬ 
tial wa ve an alysis of the yn —» rfn reaction from the Giessen 
group ifTT^ . Knowledge about the yn couplings, A "^2 
A'^ 12 ^ for the Af(1720)F’i3 is very poor, and different groups 
have given very different predictions. In the SU(6)®0(3) sym¬ 
metry limit, we predict the should be zero, which is com¬ 
patible with the analysis of the Kent group Iz^ . More studies 
are needed to clarify these puzzles about the Af(1720)Pi3. 



FIG. 13: (Color online) Total cross section as a function of the 
c.m. energy IV for the yn —> 7 r®n reaction. Data are taken 
from (3]. The results by switching off the contributions from 
M1650)5u, A( 1620 )S 3 i, and t channel; and the partial cross sec¬ 
tions for A(1232)P33, A^(1535)5 u, iV(1520)Di3, A'(1720)Pi3 and u 
channel are indicated explicitly by different legends in the figure. 


IV. SUMMARY 

In this work, we have studied neutral pion photoproduction 
on nucleons within a chiral quark model. We have achieved 
reasonable descriptions of the data from the pion production 
threshold up to the second resonance region. 

The roles of the low-lying resonances in the reactions 
were carefully analyzed. We found that; (i) A( 1232)^33, 
N{1535)S 11 , Af(1520)Di3, and Af(1720)F’i3 play crucial roles 
in both yp —» 7r“p and yn iiPn reactions. The A(1232)P33 
resonance not only plays a dominant role around the first res¬ 
onance region, but also contributes up to the third resonance 
region. Both Af(1535)5' n and Af(1520)Di3 paly crucial roles 
around the second resonance region. The second bump struc¬ 
ture around Ey - 700 MeV in the cross section receives ap¬ 
proximately equal contributions from these two resonances. 
V(1720 )Pi 3 might play a crucial role in the third resonance 
region. It might be responsible for the third bump structure 
in cross section, which should be further investigated due to 
our relatively poor descriptions of the polarization observ- 
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FIG. 14: (Color online) Cross section ratio (T,Jcrp between the reac¬ 
tions yn —> 7r°n and yp p as a function of the center-of-mass 
energy W. Data are taken from Ref. a. 


ables in this energy region, (ii) Furthermore, obvious evi¬ 
dence of Ai( 1650)5 11 and A( 1620)5 31 is also found in the re¬ 
actions. They notably affect the cross sections and the po¬ 
larization observables from the second resonance region to 
the third resonance region, (iii) The u- and f-channel back¬ 
grounds play a crucial role in the reaction as well. The ii 
channel has strong interferences with the resonances, such as 
A(1232)F33, A^(1535)5 ii and Ai(1520)Di3. By including the 
f-channel vector-meson exchange contribution, the descrip¬ 
tions of the data in the energy region Ey = 600 ~ 900 MeV 
are improved notably, (iv) No obvious evidence of the other 
resonances, e.g., N(1700)Di3, N(1675)Di5, A(1700)D33 and 
N( 1680 )Fi 5 , was found in the reactions. 

Furthermore, the helicity couplings for the main reso¬ 
nances, A(1232)F33, N(1535)5 ii, N(1520)Di3, N(1720)Pi 3, 
N(1650)5ii and A(1620)53i, were extracted from the re¬ 
actions. We found that: (i) Our extracted helicity ampli¬ 
tudes of A(1232)F’33 and N(1535)5 n are in good agreement 
with the PDG average values and the results of other groups, 
(ii) The yp coupling for N( 1650)5 n extracted us is in 
good agreement with the results from the SAID 11231 - 1^ and 
BnGa 11^ u^ . However, properties of the yn coupling for 
A(1650)5 11 are still controversial. Our analysis of the recent 
data of the yn —> uPn reaction indicates a small negative yn 
coupling for A( 1650)5 n. Its si gn i s op posite to that of other 
analyses of the yn —> pn data 1^ 1^ . (iii) We obtain 
a large helicity coupling for A( 1620)5 31 , but it is very close 
to the recent analysis from the BnGa group (iv) 

We give smaller helicity couplings for A(1520)Di3, which 
are compatible with the PDG values at the 30% level, (v) 
The helicity couplings and for A(1720)Pi3 ex¬ 
tracted by us are consi stent with the quark model predic¬ 
tions by Li an d Cl ose and the analysis of the 

Giessen group CH. We find a small positive helicity cou¬ 
pling Ay 2 for A(1720 )Pi 3, and the should be zero in the 
SU(6)®0(3) symmetry limit. 

Finally, it should be pointed out that (i) the width of 
A(1720 )Pi 3 extracted by us is notably narrower than the es- 


TABLE VI: The expressions of ^ in Eq. 
nances. Here we have defined TC = 


for various reso- 
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timated values from the PDG, however, our result is in good 
agreement with those extracted from the n^p —> K^h reaction 
by Saxon et al. ifTdl] . To confirm the properties of A(1720)Pi3, 
a study of the n^p —> K^K reaction is needed, (ii) Further¬ 
more, a more realistic correction of the FSI for neutral pion 
photoproduction on quasi-free neutrons hopefully will be ob¬ 
tained in future. Then the sign problem of the yn coupling 
A” 12 of N( 1650)5 n could be clarified in the yn n^n reac¬ 
tion, which seems to be crucial to uncover the puzzle of the 
narrow structure around W - 1.68 GeV observed in the ex¬ 
citation function of 77 production off quasi-free neutrons. If 
the yn coupling Aj^^ of the A( 1650)5 n is negative, then the 
narrow structure in the yn pn reaction would no longer be 
explained by the interference effects between A(1535)5 n and 
A(1650)5 11 . 
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FIG. 15: (Color online) Beam asymmetry of the yn —> n^n reaction as a function of scattering angle. Data are taken from ca. The first and 
second numbers in each figure correspond to the photon energy Ey (MeV) and the nN center-of-mass energy W (MeV), respectively. 


TABLE VII: Branching ratio of the resonances used in the calculation. 


Resonance 

A( 1232)^33 

A( 1620)5 31 

A(1535)5i, 

A( 1650)511 

A(1520)Di3 

A(1720)Pi3 

^7tN 

1.0 

20 ~ 30% 

35 ~ 55% 

50 ~ 90% 

55 ~ 65% 

60 ~ 90% 
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